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PREFACE 


This  Memorandum  provides  physical  data  useful  for  the  development 
of  engineering  design  criteria  for  space  vehicles  grazing  or  entering 
the  upper  atmosphere  of  Mars.  It  represents  an  extension  of  earlier 
work,  notably  P-2639,  Engineering  Model  Atmosphere  of  Mars,  to 
hl^er  altitudes  Into  the  aerodynamic  flight  regime  of  orbiting 
satellites  and  lifting  entry  vehicles.  Hie  work  was  performed  under 
Contract  AF  49(638)-700. 
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ABSTRACT 

Probable  upper  and  lower  liinlts  are  here  calculated  for  the 
distribution  of  atmospheric  pressure  and  mass  density  In  the  upper 
atmosphere  of  Mars .  The  results  extend  an  earlier  engineering  model 
atmosphere  from  an  altitude  of  150  km  to  one  In  excess  of  2500  km  above 
the  planetary  surface.  They  derive,  In  part,  from  a  recent  analysis 
by  J.  W.  Chamberlain  of  the  probable  thermal  regime  In  the  upper 
atmosphere  of  planets,  while  taking  Into  account  our  present  uncertainties 
about  the  lower  atmosphere  of  Mars. 

The  three  self-consistent  model  atmospheres  (Tentative  Maximum, 
Tentative  Minimum,  and  Tentative  Standard)  should  bracket  actual  conditions 
In  the  Martian  atmosphere  up  to  500  km  altitude  and,  above  that,  plausibly 
represent  the  extreme  range  of  probable  conditions.  At  any  specific 
altitude  level,  the  true  dally  mean  values  prevailing  during  any  season 
over  middle  and  low  latitudes  should  fall  between  these  limits.  If  no 
specific  season,  time  of  day,  or  latitude  is  specified,  the  Tentative 
Standard  Atmosphere  will  estimate  In  orders  of  magnitude  for  the  whole 
altitude  range  those  pressures  and  densities  that  will  probably  occur 
more  often  than  not.  The  parametric  limits  presented  here  are,  however ^ 
only  two  of  an  Infinity  of  choices  permitted  by  presently  available 
factual  knowledge. 
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I.  INTRODUCTION 

Ihe  AFOSR-GE  Symposium  on  Dynamics  of  Manned  Lifting  Planetary 
Entry  In  October,  1962,  demonstrated  that  the  very  high  regions  of 
planetary  atmospheres  were  Important  for  the  system  design  of  certain 
types  of  advanced  space  vehicles.  The  aerodynamlclst  and  design 
engineer  were  demanding  data  that  the  astrophysicist  and  meteorologist 
were  reluctant  to  supply  because  of  the  almost  complete  lack  of  factual 
information. Recent  policy  announcements  Indlcate^^^  that  Mars  will 
probably  be  the  first  goal  of  manned  Interplanetary  flight,  and  that 

(3\ 

serious  studies  towards  this  objective  are  In  progress .  An  attempt 
seems  warranted,  therefore,  to  provide  the  basic  physical  data  for  the 
development  of  engineering  design  criteria  for  vehicles  grazing  or 

entering  the  Martian  upper  atmosphere. 

(4) 

In  previous  studies,  the  author  has  used  whatever  factual 
information  was  available  to  derive  limiting  model  atmospheres  of  Mars 
and  to  develop  an  engineering  model  atmosphere up  to  an  altitude  of 
150  km  above  the  surface.  Wegener, for  example,  has  used  some  of 
these  scientific  data  to  describe  the  aerodynamic  flight  regimes  of 
Mars,  but  has  pointed  out  that,  even  at  altitudes  of  200  km,  the 
atmosi^iere  Is  aerodynamlcally  dense  and  that  a  fll^t  vehicle  would 
still  encounter  the  continuum  regime.  On  the  other  hand,  Vachon^^^ 
has  Illustrated  by  numerical  examples  the  well-known  fact  that 
straight  mathematical  extrapolation  of  low-altitude  data  to  altitudes 
above  a  few  hundred  kilometers  will,  with  the  assumption  of  an 
exponential  atmosphere,  lead  to  very  serious  misconceptions  about  the 
physical  state  of  any  planet's  upper  atmosphere. 

Finally,  one  must  remember  that  we  are  still  uncertain  about  the 
upper  atmosphere  of  our  own  planet.  Only  since  the  advent  of  artificial 
satellites,  have  we  begun  to  get  reliable  quantitative  Information  about 
the  large  variations  of  temperature,  pressure,  and  density  which  occur 
In  tiie  Earth's  upper  atmosphere  as  a  function  of  solar  activity, 

/g\  (9) 

time  of  day,'  '  and  season.  '  Nevertheless,  the  following  sections 
describe  a  method  for  extending  the  probable  limits  of  Mars'  basic 
ataospharlc  parasMtars  to  altitudes  In  excess  of  2500  km.  The  author, 
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knowlng  how  speculative  his  method  is,  cautions  the  reader  that  the 
tabulated  parametric  limits  may  be  seriously  in  error;  in  fact,  they 
may  be  utterly  wrong  at  altitudes  above  500  km.  Hopefully,  observations 
of  Mars  by  unmanned  space  probes  or  from  earth-based  high-altitude 
balloons  should  soon  give  us  some  answers. 
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II.  METHOD  OF  APPROACH 

An  earlier  work^^^  emphasized  the  scarcity  of  "factual,"  numerically 
reliable  data  on  the  Martian  atmosphere.  To  avoid  relying  on  subjective 
"best"  values,  that  work  used  optimization  methods  resulting  In  model 
atmospheres  that  constituted  probable-extreme  limits  for  the  basic 
atmospheric  parameters  up  to  some  150  km  In  altitude.  Actual  conditions 
over  middle  and  low  latitudes  on  Mars,  regardless  of  time  of  day  or 
season,  could  then  be  expected  to  fall  within  these  limits. 

A  number  of  studies  are  now  available  in  the  recent  literature,  e.g., 
the  excellent  work  by  Zimmerman  and  Jones and  Wegener, using 
similar  approaches  to  convert  the  astronomical  and  atmospheric  data 
for  Mars  Into  parameters  for  engineering  and  the  description  of  fllgiht 
regimes.  But,  as  mentioned,  an  extrapolation  of  such  quantities  to 
higher  altitudes  cannot  be  made  mathematically  alone;  It  must  take 
Into  account  a  wide  variety  of  complex  but  possible  physical  processes 
In  the  upper  atmosphere  of  a  planet. 

Ihls  Memorandum  will  use  extensively  the  deductive  model  of  the 
upper  atmosphere  of  Mars  developed  by  Chamberlain. His  basic 
Ideas  about  radiative  processes  In  the  Martian  thermosphere  were 
combined  with  ours  concerning  optimization;  the  results  yielded  a 
reasonable  extension  of  the  earlier  engineering  model  atmosphere^^^ 
to  altitudes  beyond  2300  km. 

In  brief,  our  reasoning  was  thus:  Chamberlain  analyzes  the  probable 
consequences  of  the  absorption  of  ionizing  and  dissociating  far-ultra¬ 
violet  solar  radiation  in  the  Martian  thermosphere  above  the  mesopause. 

He  derives  the  probable  conductive  heat  flux  through  the  upper  atmosphere 
and,  through  an  Iterative  technique,  bases  the  properties  of  the 
thermosphere  (the  region  where  the  temperature  begins  to  Increase  with 
altitude  due  to  solar  energy  absorption  through  photo- ionization  and 
photodissociation)  on  the  properties  of  the  mesopause.  In  consequence. 
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he  has  to  locate  the  mesopause  arbitrarily  at  an  altitude  where  a 
number  of  complex  conditions  may  be  fulfilled“-specifically  at  a 
density  value  where  the  conductive  heat  flux  is  completely  radiated 
away,  but  where  the  temperature  is  quite  uncertain. 

We  now  argue  that  far-ultraviolet  solar  radiation  will  penetrate 
from  above  the  atmosphere  down  to  a  certain  depth,  i.e.,  the  level  of 
the  mesopause,  independent  of  the  pressure  and  density  structure  of 
the  atmosphere  below  this  level.  The  penetration  depth  will  depend 
on  the  amount  of  atmospheric  mass  above  the  mesopause,  and  the  absorp¬ 
tion  characteristics  in  the  thermosphere,  i.e.,  the  molecular  and 

(12) 

atomic  composition.  Further,  Bates  has  shown  that  rates  of  gain 

and  loss  of  thermal  energy  in  a  planetary  thermosphere  are  rather 

insensitive  to  the  temperature  and  influenced  only  indirectly  through 

(13) 

the  atmospheric  density.  Barth  has  pointed  out  that  the  compo¬ 
sition  Of  the  upper  atmosphere  of  Mars  may  be  similar  to  that  of 
Earth  in  that  nitrogen  will  be  the  major  component  responsible  for 
radiation  effects  at  high  altitudes.  This  is  borne  out  by  Chamberlain's 
quantitative  data  (see  Ref.  11,  Fig.  2). 

To  specify  the  probable  range  of  altitudes  for  the  Martian  meso¬ 
pause,  we  use  Chamberlain's  discussion  of  thermal  heat  flux  to  deter¬ 
mine  the  amount  of  mass  penetrated  by  solar  far-ultraviolet  radiation 
down  to  the  mesopause.  The  exact  altitude  above  the  planetary  surface 
will  critically  depend  on  assisnptions  made  both  about  the  values  of 
atmospheric  parameters  near  the  surface  of  Mars  and  about  their 
vertical  distributions  in  the  troposphere  and  mesosjAiere.  But 
regardless  of  any  initial  starting  values,  the  atmospheric  pressure 
at  the  mesopause  will  be  unaffected  by  the  temperature  distribution 
above  it,  and  Chamberlain's  conditions  for  the  location  of  the  meso¬ 
pause  will  be  fulfilled  somewhere  in  the  Martian  atmosphere.  From 

(4) 

our  limiting  model  atmospheres  we  can  find  not  only  the  altitudes 
where  these  conditions  are  fulfilled  for  upper  (as  well  as  lower) 
extreme  limits,  but  also  the  temperatures  associated  with  these  levels. 

Once  the  Martian  mesopause  has  thus  been  located  for  both  upper 
and  lower  limiting  conditions,  these  limits  could  theoretically  be 
extended  upward  throu^  the  thermosphere  and  exos^ere  to  interplanetary 
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space.  But  the  energy  state  of  a  planetary  atmosphere  has  a  variety 
of  well-known  complications.  For  convenience^  studies  of  the 

Earth's  upper  atmosphere  have  divided  it  on  the  basis  of  those  compli¬ 
cations  that  are  most  significant  in  terms  of  thermal  state  and 
behavior.  Likewise^  plausible  assumptions  and  restrictions  can  be 
onade  by  so  treating  the  Martian  upper  atmosphere. 

The  Appendix  gives  a  more  detailed  discussion  of  the  applied 
concepts  and  mathematical  details.  From  various  assumptions^  we  were 
able  to  extend  the  altitude  variation  of  atmospheric  parameters  up  to 
2600  km  above  the  planetary  surface,  getting  both  upper  and  lower 
probable  limits  even  at  these  high  altitudes. 


III.  PARAMETRIC  RESULTS 


Numerically  applying  the  concepts  described  In  the  Appendix  produced 
a  series  of  model  atmospheres  for  Mars,  each  self-consistent  and  physically 
valid  for  a  particular  configuration  of  planetary  and  atmospheric  data. 
Since  the  objective  was  to  derive  upper-  and  lower-probable  limits  for 
the  variation  of  pressure  and  density  in  the  upper  atmosphere  of  Mars, 
three  models  that  reasonably  fulfill  this  requirement  were  selected 
and  labelled  as  follows: 

TENTATIVE  MAXIMUM  —  probable  upper  limit, 

TENTATIVE  STANDARD  —  probable  average, 

TENTATIVE  MINIMUM  —  probable  lower  limit . 

The  principal  characteristics  of  these  model  atmospheres  are 
illustrated  in  Figs.  1,  2  and  3,  and  are  summarized  below.  Values  of 
atmospheric  pressure  and  density  from  the  planetary  surface  to  2800  km 
altitude  are  tabulated  in  Table  1.  A  more  detailed  discussion  is  found 
in  Appendix  D,  and  Table  6  gives  the  construction  parameters  while 
Table  7  gives  the  thermal  regimes  with  an  example  of  a  representative 
molecular  mass  distribution. 

TENTATIVE  MAXIMUM 

This  model  atmosphere  represents  extreme  conditions  that  may  be 
encountered  if  the  lower  portion  of  a  dense  Martian  atmosphere  is  heated 
with  maximum  efficiency  by  absorption  of  solar  radiation  in  both  the 
troposphere  and  the  mesosphere.  The  mesopause  level  is  found  at  an 
altitude  of  about  215  km  and  at  a  temperature  of  some  At  higher 

altitudes,  the  atmosphere  is  again  absorbing  solar  energy  with  high 
efficiency,  causing  dissociation  of  the  constituent  gases.  It  is 
assumed  —  for  this  extreme  case  —  that  the  molecular  scale  temperature 
could  reach  4400°K  at  an  altitude  of  2700  km.  This  high  molecular 
scale  temperature  would  correspond  to  a  kinetic  temperature  ranging  from 
1100°  with  a  mean  molecular  mass  of  7  (e.g>;  ionization  as  well  as  complete 
dissociation,  and  helium,  hydrogen,  or  oxygen  the  dominant  species) 
to  4400°K  (without  any  dissociation). 
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This  model  atmosphere,  in  terms  of  the  altitude  variation  of 
pressure  and  mass  density,  might  be  found  on  a  Martian  summer  day  in 
equatorial  latitudes  during  perihelion.  The  critical  level  of  escape 
is  to  be  found  at  extremely  hi^  altitudes,  in  excess  of  2500  km. 

Thus,  the  noticeable  atmosphere  reaches  much  farther  from  the  planetary 
surface  than  it  does  on  Earth. 

TENTATIVE  MINIMUM 

This  model  atmosphere  is  valid  for  extreme  conditions  of  a  cold, 

thin  Martian  atmosphere .  The  mesopause  level  is  located  near  91  km 

altitude  at  a  temperature  of  about  101°K.  A  gradual  rise  in  molecular 

scale-temperature  leads  to  values  of  906°K  near  500  km,  and  1600°K 

near  1500  km.  The  escape  level  is  between  590  km  and  635  km  at  a 

molecular- scale  temperature  of  985°K  to  1020°K.  This  corresponds 

to  an  approximate  kinetic  temperature  of  900°K  with  dissociation  to  a 

o 

molecular  weight  of  26  or  27,  or  to  one  of  1000  K  without  any  appreciable 
dissociation.  It  might  be  found  on  a  Martian  winter  night  in  middle 
latitudes  near  aphelion. 

TENTATIVE  STANDARD 

This  model  atmosphere  was  selected  as  being  representative  of 
M’^rtian  conditions  thought  to  average  by  most  authors.  The  mesopause 
occurs  at  an  altitude  of  144  km  with  a  temperature  of  162. 5°K.  The 
molecular  scale  temperature  rises  gradually  to  1990°K  near  1500  km 
and  continues  to  increase  to  almost  2300°K  at  approximately  2700  km. 

From  above  1500  km  the  true  kinetic  temperature  probably  remains 
almost  constant  near  1100*^K,  the  molecular  scale-temperature  rising 
because  of  molecular  dissociation,  with  consequent  decrease  of  the 
mean  molecular  mass  to  values  of  14  near  2600  km. 

Depending  on  the  precise  variation  of  gravitational  acceleration 
with  altitude,  primarily  the  latitude  concerned,  the  critical  level  of 
escape  is  between  1540  km  and  1760  km.  The  molecular  scale -temperature 
there  is  between  200o\  and  2070°K,  e.g. ,  corresponding  to  a  mean 
molecular  wel^t  of  15.7  for  a  kinetic  temperature  of  1100°K. 
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This  Tentative  Standard  Model  Atmosphere  probably  represents 
no  actual  conditions  found  on  Mars  at  any  one  time,  place,  or  altitude. 
But  it  should  indicate  the  approximate  values  of  pressure  and  temperature 
likely  under  conditions  averaged  over  day,  season,  latitude,  and  orbit 
position.  As  with  model  atmospheres  for  Earth, its  primary  use 
will  be  for  preliminary  engineering  designs. 

GENERAL  APPLICATION 

These  tentative  model  atmospheres  will  be  useful  if  their  purposes 
are  kept  in  mind.  The  Tentative  Maximum  and  Tentative  Minimum  models 
will  give  valid  numerical  limits  for  values  of  atmospheric  pressure  and 
mass  density  at  altitudes  ranging  from  0  km  to  some  2500  km.  At  any 
specific  altitude,  the  true  values  prevailing  at  any  time  or  season 
over  middle  and  low  latitudes  should  fall  between  these  limits.  The 
Tentative  Standard  Atmosphere,  on  the  other  hand,  will  give  those 
order-of-magnitude  estimates  of  pressure  and  density  for  the  whole 
altitude  range  —  estimates  that  will  be  close  to  reality  more  often 
than  not,  if  no  specific  season,  time  of  day,  or  latitude  is  specified. 
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IV.  DISCUSSION 

Hie  numerical  data  can  only  be  as  reliable  as  the  many  assumptions 
that  preceded  the  computation.  Unfortunately,  there  Is  not  enough 
factual  Information  available  to  assess  In  any  way  their  reliability  or 
probability.  Nevertheless,  the  method  used  here  will  probably  be  less 
In  error  than  any  straight  extrapolation  of  lower-altitude  values  from 
an  Isothermal  upper  atmosphere  without  dissociation,  l.e.,  based  on  an 
exponential  variation  of  density,  with  altitude  above  the  mesopause.  This 
Is  so,  however,  only  because  the  assumption  of  a  constant- temperature/ 
molecular-mass  regime  Is  untenable  from  the  point  of  view  of  planetary 
physics,  and  would  lead  to  the  maximum  possible  errors  In  hlgh-altltude 
densities . 

In  essence,  this  analysis  extends  that  of  Chamber laln^^^^  to  allow 
for  our  uncertainties  concerning  the  lower  atmosphere  of  Mars.  In 
addition,  little  Is  said  here  about  the  molecular  and  atomic  constituents 
of  the  Martian  thermosphere  and  exosphere  above  some  500  km.  Earlier 
work  by  Yanow^^^^  Indicated  how  complex  are  the  processes  which  may 
occur  I'*  this  region  and  lead  to  Ionized  layers  analogous  to  the  Earth's 
Ionosphere.  Rasool^^®^  recently  criticized  Chamberlain's  analysis  of 
the  Martian  Ionosphere  for  having  an  escape- level  kinetic  temperature 
so  low  (1100°K)  that  it  could  not  allow  for  diurnal  variations  nor 
for  Increases  associated  with  solar  activity;  hence,  its  electron-density 
distribution  might  be  an  upper  limit.  Rasool  does  not  offer  an  alternative 
approach,  however,  beyond  pointing  out  that  temperature  variations  in  the 
Earth's  upper  atmosphere  are  known  to  be  possibly  as  high  at  500°K  as 
a  function  of  local  time  and  solar  activity.  Nor  does  Rasool  seem  to 

(19) 

have  considered  the  role  of  diffusive  equilibrium  at  these  altitudes, 
a  role  that  tends  to  fortify  Chamberlain's  numerical  results. 

Further  critical  analysis  of  the  exosphere  raises  questions  as 
to  how  far  from  the  surface  does  the  atmosphere  continue  to  rotate  with 
the  planet  as  a  solid. Hiough  the  answer  to  this  question  depends 
upon  the  magnitude  of  a  Martian  magnetic  field  and  upon  its  interaction 
with  solar  wind,  the  present  study  Implicitly  assumed  that  the  Martian 
atmosphere  does  rotate  with  the  planet  up  to  the  altitudes  considered. 
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V.  CONCLUSIONS 


An  infinity  of  self-consistent  model  atmospheres  of  Mars  could 
be  constructed  from  presently  available  factual  knowledge.  Though  the 
parametric  limits  presented  here  are  just  two  from  this  infinity,  they 
are  probably  rather  reliable  for  altitudes  from  the  surface  up  to  some 
500  km.  Between  500  km  and  1500  km,  they  plausibly  represent  the 
range  of  probable  conditions.  For  still  higher  altitudes,  the  upward 
extension  at  a  constant  exospheric  temperature  (kinetic)  can  presently 
be  considered  only  as  a  best  guess.  This  does  imply,  in  essence,  that 
the  Martian  exosphere  is  at  an  altitude  of  1500  km  or  more. 

The  tentative  model  atmospheres  show  that  despite  the  planet's 
smaller  volume  and  mass,  the  sensible  atmosphere  extends  to  propor¬ 
tionally  much  higher  altitudes  than  on  Earth.  This  anomaly  is  caused 
by  the  atmospheric  dissipation  rate's  depending  not  on  the  gravita¬ 
tional  acceleration,  but  on  the  gravitational  potential;  it  is  the 
latter  that  determines  the  critical  escape-velocity  of  the  atoms  and 
molecules,  as  well  as  the  altitude  variations  of  pressure  and  density. 
This  must  always  be  remembered  when  trying  to  "compare"  the  physical 
characteristics  of  different  planetary  atmospheres  at  various 
geometric  altitudes  above  the  surfaces. 

Numerical  comparison  with  Earth, for  example,  immediately 

demonstrates  that  artificial  satellites  will  experience  considerably 

higher  drag  effects  while  orbiting  Mars  than  at  similar  altitudes 

above  the  Earth's  surface.  To  show  this.  Figs.  1,  2,  and  3  contain 

an  upper  scale  of  representative  Earth  altitudes  for  each  parameter. 

These  Earth  scales  are  based  on  the  U.  S.  1962  Standard  Atmosphere 

below  400  km.  Above  this  altitude,  they  are  based  on  satellite  drag 

data  and  reflect  the  variations  of  temperature,  pressure,  and  density 

as  a  function  of  local  time  and  solar  activity  in  the  Earth's  upper 

atmosphere.  In  tetms  of  the  planetary  gravitational  fields,  on  the 

other  hand,  a  geometric  altitude  of  2700  km  on  Mars  corresponds  to 

one  of  about  630  km  on  Earth;  and  we  find,  for  example,  that  atmos- 

-12 

pheric  pressures  on  both  planets  have  dropped  to  about  10  of 
their  respective  surface  pressures  at  this  radial  distance.  The 
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density  scale-height  is  a  governing  parameter  for  planetary  entry. 

For  meaningful  comparison  with  Earthy  however^  we  must  consider  it  in 
terms  not  of  geometric-altitude  intervals  but  of  the  relative  decrease 
of  atmospheric  density  over  corresponding  differences  of  gravitational 
potential  surfaces. 

Finally,  as  with  any  standard  atmosphere, it  is  well  to 
remember  that  the  values  falling  between  the  envelope  curves  of  any 
of  the  parameters  could  not  possibly  be  encountered  at  all  altitudes 
at  any  given  location  and  time.  In  other  words,  the  maximum  spread 
is  not  representative  of  actual  variability  at  any  altitude,  but 
merely  provides  for  extreme  conditions  that  could  exist  somewhere 
over  nearly  the  entire  globe  of  Mars.  At  no  time  is  it  likely  that 
the  atmosphere  will  assume  the  lower  (or  upper)  limit  distribution 
throughout  its  depth,  or  through  any  major  segment  thereof. 

On  Earth,  for  instance,  the  warmest  surface  layers  in  the 
tropics  are  known  to  be  associated  with  the  coldest  (and  highest) 
tropopause  levels.  In  Fig.  1,  the  temperature  envelopes  indicate 
the  possibility  of  an  analogous  phenomenon  —  the  Martian  mesopause  — 
governing  the  vertical  distribution  of  atmospheric  density.  But  any 
specific  altitude  distribution  of  atmospheric  parameters  throughout 
the  atmosphere  within  the  permissible  reuiges  will  also  evidence  this 
characteristic.  Hence  we  have  difficulty,  without  better  observa¬ 
tional  data,  in  providing  physically  significant  mean  or  standard 
altitude  distributions,  or  in  narrowing  the  parametric  limits. 

In  fact,  very  recent  spectroscopic  observations  have  led 

(22) 

Kaplan'  to  conclude  tentatively  that  the  surface  pressure  on  Mars 
may  be  as  low  as  20  mb.  Such  an  atmosphere,  probably  consisting 
mostly  of  CO^  and  argon,  would  fall  outside  of  even  the  lower 
probable  limits  derived  here.  This  is  a  further  illustration,  if  we 
needed  one,  that  present  factual  knowledge  does  not  permit  us  to 
specify  any  degrees  of  probability  to  the  present  work. 
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Appendtx  A 

MATHEMATICAL  TECHNIQUE 


A  program  has  been  developed  for  the  IBM  7090  Computing  Machine, 
which  rigorously  evaluates  the  hydrostatic  equation  under  the  assump¬ 
tion  that  the  atmosphere  behaves  as  an  ideal  gas.  While  this  proce¬ 
dure  would  be  trivial  for  an  exponentially  decreasing  atmosphere  in 
an  inverse-square,  central  gravitational  field,  an  actual  planetary 
atmosphere  such  as  Mars '  requires  a  solution  that  takes  into  account 
the  empirical  data  we  have  as  well  as  that  which  we  lack.  This  we 
can  achieve  by  so  designing  the  computing  program  as  to  permit  substi 
tution  of  known,  composite  variables  for  unknown  data  at  any  altitude 
level  with  explicit  requirements  for  necessary  assumptions. 

Let  the  perfect  gas  law  be: 

* 

p  R  T 

p  -  -  -  n  k  T  ,  (1) 


where 


m 

_2 

p  ■  atmospheric  pressure  (in  dynes  cm  ), 

_3 

p  •  mass  density  (in  g  cm  ), 

T  ■  true  kinetic  temperature  (in  °K), 
m  ■  mean  molecular  mass, 

3 

n  ■  number  density  of  particles  (per  cm  ) 
* 

K  "  universal  gas  constant,  and 
k  ■  Boltzmann  constant. 


where  the  hydrostatic  equation  to  be  integrated  is 

‘‘f  •  •  *(r)  "(r) 

where 

g^^^  ■*  gravitational  acceleration  at  altitude  z 
(in  cm  sec  and 

z  ■  geometric  altitude  above  MSL  or  mean  planetary 
surface  (in  cm). 
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The  Integration  to  obtain  p  and  p  at  any  altitude' z  can  be 
performed  through  substitution  for  the  variables  z,  T,  and  m  as 
detailed  below. 


POTENTIAL  HElfflT 

A  potential  height  h  is  defined  as 

z 

.  ’O  Z  I* 

h  -  —  *  —  r  g 


8. 


1 

-  ! 
%  ° 


(2) 


dz, 


where 


h  *  potential  height  above  planetary  surface  (in  cm), 

i  s  gravitational  potential  at  altitude  z 
2  -12-2 

(in  erg  gm  =  cm  sec  ) ,  and 

g^  ■  gravitational  acceleration  at  zero  altitude 
(constant  for  a  specific  latitude). 

If  empirical  data  on  the  distribution  of  g^  with  planetary 
latitude  on  a  non-spherical,  rotating  planet  are  available,  the  true 
planetary  radius  at  a  given  latitude  can  be  replaced  by  an  "effective 
planetary  radius,"  so  that 


g 


(2) 


g. 


°  (Ro  ^ 


(4) 


where  ■  effective  planetary  radius  at  a  specific  latitude;  thus, 

g  as  well  as  g,  .  are  valid  at  this  chosen  latitude.  This  is 
*o  »(z)  ^21) 

successful  on  Earth  and  yields,  for  the  gravitational  accelera¬ 
tion  at  great  altitudes  over  specific  latitudes,  numerical  values 
that  are  as  accurate  as  approximations  based  on  the  expansion  of 
the  gravitational  potential  in  infinite  series  of  spherical  harmonics 
with  empirical  values  of  the  second  harmonic  coefficient. 

Lacking  empirical  data  for  the  latitude  dependence  of  g^  on  Mars, 
we  chose  the  following  planetary  radii  to  reflect  upper  and  lower 
probable  limits  and  a  probable  standard  condition: 


Table  2 


PRINCIPLE  PLANETARY  PARAMETERS 


Sym¬ 
bol  Parameter 

Units 

- 1 

Lower  Limit 

Tentative 

Standard 

Upper  Limit 

Effective 

R  planetary 

cm 

3440  X  10^ 

3375  X  10^ 

3322  X  10^ 

radius 

Effective 
sur face 
acceleration 

-2 

cm  sec 

360 

375 

390 

of  gravity 

Representa- 
M  tive  plane- 

gm 

26 

6.389  X  10 

6.406  X  10^^ 

6.455  X  10^^ 

tary  mass 

Note  that  the  values  for  planetary  masses  are  illustrative  only^  and 
would  not  necessarily  represent  the  non-rotating  inertial  mass. 

Using  "effective  planetary  radius"  to  account  for  flattening 
and  rotation  of  Mars^  we  can  integrate  Eq.  (3)  to  yield 


R  h 
o 

R  -h 
o 


or 


R  z 
0 

R  +z 
o 


(5) 


The  resulting  potential  height  h  (in  cm)  will  be  an  artificial  measure 
(particular  to  Mars)  of  the  increase  of  gravitational  potential  and 
the  decrease  of  gravitational  acceleration  with  altitude  for  a  specific 
latitude. 


APPARENT  GRAVITATIONAL  FIELD 


If  the  physical  characteristics  of  a  planet  are  known  with  higher 
accuracy^  the  concept  of  potential  height  can  be  further  expanded  to 
represent  the  forces  acting  at  different  latitudes  upon  an  atmosphere 
rotating  with  the  planet.  In  Eq.  (3)^  the  approximate  exfhressions  for 
the  gravitational  potential  and  its  gradient  become: 


CM 


(2 .9) 


1  + 


2 


(1-3  sln^cp) 


n‘  2  2 

+  —  r  cos  cp  , 

2 


(6) 
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where 


and 


r 


(1  -  f  sin^  (p)  +  z  , 


(7) 

(8) 


where 

r  *  radius  vector, 

G  ■  universal  gravitational  constant, 

M  ■  planetary  mass, 

=  equatorial  radius, 

J2  ■  second  harmonic  coefficient, 
n  “  angular  rotation, 
cp  “  planetary  latitude,  and 
f  •  flattening. 

It  is  evident  that  with  these  input  quantities,  the  computer 
program  becomes  generally  applicable  to  any  planet.  This  routine 
was  used  to  compare  results  at  different  latitudes  at  higher  altitudes 
(see  Appendix  F).  Within  the  still  existing  uncertainties  of  the 
physical  characteristics  of  Mars,  both  routines  are  obviously  more 
accurate  than  input  data  warrant. 


MOLECULAR  SCALE  TEMPERATURE 


Since  the  mean  molecular  weight  is  also  a  function  of  altitude, 
a  second  substitution  is  necessary,  defining  a  "Molecular  Scale 
Temperature"  T^  as 


where 


(9) 


m  ■  mean  molecular  mass  at  zero  altitude  ■  constant, 
o 

m  ■  f(z)  or  f(h), 

T  •  f(z)  or  f(h),  and 
m 

T  ■  f(z)  or  f(h). 

With  the  substitutions  of  molecular  scale  temperature  and 
potential  height,  Eq.  (2)  becomes 


dp  “  -  gjj  P  dh  . 


(10) 
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With  the  help  of  Eqs .  (1)  and  (9) ,  it  can  be  Integrated  with  certain 
restrictions  on  the  variations  of  molecular  mass  and  kinetic  temperature 
with  altitude,  yielding 


P(h2)  ^’(hl) 


g  m 
o  o 


R 


dh 

t” 

m 


(11) 


Specifically,  if  the  variation  of  molecular  scale  temperature  with 
potential  height  is  either  zero  or  constant  in  a  certain  height  interval, 
that  is,  where 


d  T  -  Ldh 
m 


(12) 


but  where 


ST 

m  , . .  m 


L  =  —  with  - =  0, 

dh  dh 


then  (for  L  =  constant); 


"(hp  "(hp  ^ 


m(h^) 


m  g 
o**o 

R*L 


for  L  ■  0  (i.e.,  “  constant  between  hj^  and  h^); 

m  g 
■  R*T 


^’(h^)  “  ‘’(hp  ® 


(13) 


(14) 


It  is  now  possible  to  evaluate  Eq.  (1)  at  all  levels  of  h  and  z, 
and  to  derive  mass  densities  p  for  given  values  of  the  ratio  T/m.  In 
addition,  any  combination  of  known  and  unknown  parameters  in  Eq.  (9) 
will  specify  what  assumptions  must  be  made  about  the  missing  information. 

The  computing  program  allows  flexibility  with  respect  to  acceptance 
of  known  parameters  and  conditions  inherent  in  Eqs.  (1),  (3),  (5),  (9), 
and  (11)  as  limited  by  conditions  in  Eq.  (12)  ,  and  can  be  called  upon  to 
provide  output  of  additional  quantities,  such  as  scale  heights,  columnar 
mass,  number  densit’  and  other  auxiliary  quantities.  As  described  in  the 


following  sections,  uncertainties  in  a  variety  of  parameters  for 
Mars  necessitated  the  employment  of  maximization  and  iteration 
techniques  to  develop  a  series  of  self-consistent  model  atmospheres, 
representative  of  limiting  conditions  in  various  altitude  Intervals. 
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Append  lx  B 
LOWER  ATMOSPHERE 

(4) 

An  earlier  study  calculated  realistic  upper  and  lower  limits 
for  the  permissible  ranges  of  temperature,  pressure,  and  density  of 
the  Martian  atmosphere  at  altitudes  below  the  then-unknown  mesopause. 
It  documented  In  detail  published  justifications  for  the  required 
choices  among  available  physical  data.  The  present  work  extends 
these  limits  upwards ,  following  (with  one  exception)  the  procedures 
of  another,  earlier  study. 

The  exception  Is  a  result  of  uncertain  values  of  Martian  mass, 
diameter,  and  flattening,  leading  to  further  uncertainties  about  the 
precise  variation  of  gravitational  acceleration  with  altitude.  For 
variation  In  the  gravitational  potential  cannot  be  neglected  at 
altitudes  above  some  150  km  without  Introducing  appreciable  errors. 

In  order  to  achieve  the  necessary  precision,  the  numerical  Inte¬ 
grations  were  performed  In  terms  of  "potential  height"  and  effective 
planetary  radius,  or  Indirectly  In  terms  of  apparent  gravitational 
potential,  as  the  previous  section  explained.  But  for  this  one 
exception,  the  tentative  model  atmospheres  correspond  below  150  km 

to  limiting  envelopes  published  earlier,  and  the  same  meteorological 

(4) 

slt\iations  apply  as  discussed  there. 

Table  3 

BASIS  OF  LOWER  ATMOSPHERIC  DATA 


Present  Work 

Earlier  Work  (with  Reference) 

Tentative  Minimum 

Model  II  Lower  Limit qv 
Lower  Limit  Model  (Engineering  Atmosphere)'  ' 

Tentative  Standard 

Model  II  Mean^^^ 

Tentative  Maximum 

Upper  Limit  Model  (Engineering  Atmosphere) 
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Append  lx  C 
MESOPAUSE 


Chamberlain^  locates  the  Martian  mesopause  where  the 

2 

conductive  heat  flux  is  about  0.07  ergs /cm  sec,  the  number  density 

"3 

is  1.2  X  10  cm  ,  and  the  temperature  approximately  76°K.  Based 

(23) 

on  the  model  atmosphere  of  Goody  for  heights  up  to  100  km,  and 
an  adiabatic  temperature  decrease  above  114  km,  he  obtains  a  meso¬ 
pause  altitude  of  about  130  km  for  these  conditions . 

From  Chamberlain's  data  of  mesopause  temperature  and  number 

density,  we  can  determine  the  amount  of  coliimnar  mass  penetrated 

-3  2  -3 

down  to  this  level  as  between  3.23  x  10  gm/cm  and  3.77  x  10 
2 

gm/cm  ,  depending  on  assumptions  about  the  thermospheric  molecular 
composition  and  the  decrease  of  gravitational  acceleration  with 
altitude.  This  condition  corresponds  to  a  pressure  level  of  about 
1.259  X  lO"^  mb. 

We  now  find  the  altitude  where  this  condition  is  valid  for  the 
upper  and  lower  limits  of  our  engineering  model  atmosphere, by 
extrapolating  its  upper  limit  from  150  km  for  an  adiabatic  distri¬ 
bution,  and  the  lower  limit  for  an  isothermal  temperature  structure, 
respectively. 

As  expected,  the  mesopause  lies  at  a  rather  low  altitude  for 
the  lower  limit  of  the  engineering  model  atmosphere,  and  a  very  high 
altitude  for  the  upper  limit,  thus  bracketing  Chamberlain's  deter¬ 
mination.  However,  it  can  be  noted  that  the  temperature  is  almost 
the  same  in  both  limiting  cases,  as  required  by  the  physics  of  the 
situation.  This  result  allows  us  to  reduce  the  probable  limits  for 
extreme  conditions  in  the  upper  atmosphere  of  Mars  considerably 
beyond  any  straight  extrapolation  of  the  engineering  model  atmos¬ 
phere. 

In  addition,  we  can  determine  the  extreme  altitude  ranges  for 
certain  characteristic  processes.  Of  special  interest  to  an  under¬ 
standing  of  the  physical  state  of  a  planetary  atmosphere  are  the 

levels  of  vibrational  and  rotational  relaxation,  as  defined  by 

(24) 

Curtis  and  Goody.  Similarly,  we  can  now  estimate  the 
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approxlmate  altitudes  vdiere  solar  radiation  pressure  equals  atmos¬ 
pheric  drag  on  a  satellite  of  representative  dimension,  or  where  the 
continuum  aerodynamic  regime  changes  to  free-molecular  flow.  These 
levels  are  summarized  in  Table  4. 

(4) 

It  was  pointed  out  originally  that  the  lower  limit  of  our 
model  atmosphere  was  relevant  for  a  very  cold  region  on  Mars  with 
low  surface  temperatures,  a  cold  troposphere  and  mesosphere,  and  the 
smallest  probable  total  of  air  mass.  In  such  an  atmosphere,  the 
mesopause  will  be  as  low  as  89  km,  and  the  temperature  will  begin  to 
increase  upwards  from  there  on. 

In  contrast,  the  upper  limit  pertained  to  a  planet  or  area 
with  high  surface  temperatures  and  an  efficiently  warmed  mesosphere. 
As  a  result,  however,  the  cold  mesopause  level  will  be  as  high  as 
202  km  and  heating  due  to  absorption  of  far-ultraviolet  solar 
radiation  will  take  place  only  at  higher  altitudes. 

In  order  to  provide  also  a  tentative  standard  atmosphere,  the 
Model  II  Mean  of  Ref.  4  was  extended  upwards  isothermally,  resulting 
in  an  Intermediate  mesopause  altitude  of  144  km. 
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Appendlx  D 
THERMOSPHERE 


Although  we  have  now  arrived  at  extreme  probable  limits  for  the 
altitude  and  temperature  of  the  mesopause  on  Mars,  we  still  need  to 
know  the  precise  absorption  characteristics  of  the  thermosphere  to 
determine  the  rate  of  temperature  increase  with  altitude  above  this 
level.  As  Chamberlain^^^^  points  out,  the  initial  temperature- 
altitude  gradient  could  be  as  high  as  20° /km,  if  the  upper  atmos¬ 
phere  of  Mars  had  the  same  chemical  composition  as  that  of  Earth. 

He  finds,  though,  that  CO  in  the  Martian  upper  atmosphere  will  act 
as  an  effective  thermostat,  reducing  this  gradient  substantially. 

His  resultant  temperature  distribution  is  achieved  through  an 

iterative  process;  we  shall  assume  that  his  assumptions  are  valid. 

(191 

Nicolet,  in  his  studies  of  the  upper  atmosphere  of  Earth, 

has  shown  that  the  vertical  distribution  of  density  at  very  high 

altitudes  (above  the  thermopause)  depends  critically  on  the 

temperature  at  the  mesopause  and  above.  But  with  diffusion  playing 

the  dominant  role,  the  variation  of  the  scale  height  gradient  in  the 

thermosphere  is  due  essentially  to  a  decrease  in  molecular  mass, 

rather  than  to  an  increase  of  the  kinetic  temperature. 

He  are  then  faced  with  the  problem  of  locating  the  thermopause 

and  the  critical  level  of  escape.  This  escape  level  will  be  reached 

at  an  altitude  where  the  mean  free  path  becomes  equal  to  the  pressure 

scale  height.  The  mean  free  path  Itself  is  a  function  of  the 

particle  density  and  the  effective  collision  diameter  (p)  of  the 

particles  involved.  Where  dissociation  takes  place,  the  mean 

molecular  mass  and  the  mean  effective  collision  diameter  will  tend 

-8 

to  decrease  from  the  value  of  about  3.6  x  10  cm  for  a  mixture  of 
nitrogen,  oxygen,  and  carbon  dioxide  molecules.  Ionization  will 

tend  to  increase  the  effective  collision  diameter,  however,  and  we 

-8 

are  probably  Justified  in  retaining  a  value  of  a  ■  3.5  x  10  cm 
throughout  the  Martian  thermosphere. 
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Chamberlain  gives  a  value  of  n  ■  1  x  10^  and  T  =*  1100°K  for  the 

level  of  escape  at  an  altitude  of  1500  km;  this  corresponds  to  a 

-9 

pressure  of  about  1.5  x  10  mb. 

From  the  relation 


* 

\ 


1 


J2  TT 


* 

n 


* 

P 

*  * 

P  g 


where 


* 

\ 


mean  free  path  at  level  of  escape, 

mean  effective  collision  diameter  of  air  molecules, 
particle  density  at  level  of  escape, 
pressure  scale  height  at  level  of  escape, 
mass  density  at  level  of  escape, 
pressure  at  level  of  escape, 

gravitational  acceleration  at  level  of  escape, 
mean  molecular  mass  at  level  of  escape, 
molecular  scale  temperature  at  level  of  escape,  and 
true  kinetic  temperature  at  level  of  escape, 


it  follows  that  the  atmospheric  pressure  p  at  .the  escape  level  will 
be  given  by: 

*  * 

*  ®  ,  2 
p  =*  ■  X  3.7316  X  10  [dynes/cm  ], 

a 

•k 

and  the  atmospheric  temperature  T  by: 


*  * 

ra  g 

X  2.7035  X  10"^  [®K]. 

n  a 

Depending  on  assumptions  about  the  value  of  the  gravitational 
acceleration  at  the  surface  and  the  planetary  radius,  Chamberlain's 
escape-level  values  apparently  correspond  to 

2.2  X  10^^  <  -^  <  2.3  X  10^^  [cm"^], 


-31- 


•8 

or  effective  collision  diameters  between  2.6  x  10  cm  and 
“8 

2.8  X  10  cm  for  mean  molecular  masses  of  16  to  17. 

In  studies  of  the  upper  atmosphere  of  Earth,  satellite-drag 
observations  lead  to  values  of  mass  density  p,  and  rocket-absorption 
measurements  to  the  knowledge  of  number  density  n.  But  these  can 
be  interpreted  only  in  terms  of  pressure  and  kinetic  temperature, 
if  the  mean  molecular  mass  is  known.  Similarly,  even  if  the  thermal 
heat  flux  through  the  Martian  exosphere  is  assumed  to  be  known,  as 
determined  by  Chamberlain,  a  singular  solution  in  terms  of  a 
consistent  model  atmosphere  will  still  be  incomplete,  as  is  discussed 
in  Appendix  A. 

For  our  purpose  of  arriving  at  reasonable  probable  limits  of 
pressure  and  density  in  the  upper  atmosphere  of  Mars,  a  series  of 
model  atmospheres  was  therefore  computed  from  different  basic 
assumptions.  This  flexible  approach  allowed  for  our  ignorance  of 
the  precise  variation  with  altitude  of  the  mean  molecular  mass  for 
the  physical  conditions  inherent  in  different  starting  conditions 
near  the  planetaiy  surface. 

The  problem  of  integration  can  be  summarized  as  follows: 

We  assume  that  atmospheric  conditions  are  known  and  specified 
at  the  level  of  the  mesopause.  We  further  assume  that  either  the 
true  kinetic  temperature  in  the  thermosphere  above  the  level  of 
escape  is  known,  or  that  the  molecular  scale  temperature  can  be 
estimated.  As  Chamberlain  discusses  in  detail, however,  only  an 
iterative  process  will  yield  the  altitude  distribution  of  the 
specific  molecular  and  atomic  constituents,  and  this  only  for  an 
assumed  relative  preponderance  of  constituents,  with  assumptions 
about  their  absorbing  efficiencies  and  the  resultant  heat  flux.  On 
the  other  hand,  since  we  are  primarily  interested  in  the  distribution 
of  mass  density  and  pressure  with  altitude,  it  is  feasible  to 
integrate  the  atmospheric  equations  in  terms  of  molecular  scale 
temperature,  regardless  of  the  true  kinetic  temperature  at  levels 
between. 

Conversely,  an  assumption  must  be  made  about  the  most  likely 
distribution  of  mean  molecular  mass  m  with  altitude.  But  the 
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posslble  errors  thus  introduced  in  terms  o£  pressures  and  mass 
densities  are  of  the  order  of  factors  of  2  or  possibly  of  3  rather 
than  orders  of  magnitude.  Figure  4  shows  this  schematically. 


Fig.  4  —  Schematic  sketch  of  variation 
of  mean  molecular  mass  with  altitude 


At  the  mesopause  level,  m  will  be  near  30  (Point  A).  Some¬ 
where  high  in  the  thermosphere  where,  for  all  practical  purposes, 
outer  space  is  approached,  m  mi^t  be  as  low  as  4  (Point  B) .  In 
such  a  highly  ionized  fringe  region,  m  is  complicated  by  an 
orbiting  component  of  particles,  electrons,  and  the  flux  of 
escaping  lifter  elements.  The  true  altitude  distribution  of  m, 
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governed  by  diffusive  separation  in  the  gravitational  field,  may 
lie  somewhere  between  the  dashed  lines.  But  any  reasonable  inte¬ 
gration  path  from  A  to  B  for  a  specific  energy  flux,  i.e.,  tempera¬ 
ture  equilibrium,  will  be  unlikely  to  deviate  from  the  true  molecular 
mass  at  any  level  by  more  than  a  factor  of  3. 

From  a  set  of  plausible  model  atmospheres  over  various  inte¬ 
gration  paths,  we  finally  selected  three  that  represented  the 
limiting  probable  conditions  in  terms  of  density  and  pressure 
distribution.  In  the  next  section.  Table  6  lists  the  precise 
construction  parameters  calculated  from  the  following  primary 
considerations. 

TENTATIVE  MINIMUM 

It  was  assumed  that,  for  this  extreme  case,  minimum  heating  and 
dissociation  takes  place  in  the  thermosphere  and  that  in  the 
exosphere  the  conditions  were  equivalent  to  a  kinetic  temperature  of 
906°K,  as  suggested  recently  by  Walker  and  Jastrow,^^^^  combined  with 
a  molecular  mass  of  17.  The  slow  decrease  from  a  molecular  mass  of 
30  (below  the  escape  level)  to  17  (near  2600  km  altitude)  was  assumed 
to  be  due  to  dissociation  of  both  oxygen,  nitrogen,  and  C02» 

TENTATIVE  MAXIMUM 

Maximum  heating  was  assumed  to  occur  in  the  thermosphere, 
reaching  exponentially  a  kinetic  temperature  of  llOO^K  near  500  km 
altitude  with  a  molecular  mass  of  22,  corresponding  to  an  even 
mixture  of  nitrogen  molecules  and  oxygen  atoms.  In  the  exosphere, 
at  a  potential  height  of  1500  km  (equivalent  to  a  geometric 
altitude  of  2735  km),  it  was  assumed  that  the  molecular  scale- 
temperature  reached  4400°.  If  the  kinetic  temperature  remained 
constant  at  1100°  K  up  to  these  altitudes,  there  would  be  a  molecular 
mass  of  7  *—  a  mixture  primarily  of  atomic  nitrogen  and  a  flux  of 
atomic  oxygen  and  helium  escaping  through  the  thermosphere  and 
exosphere . 
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As  was  explained  earlier,  however,  this  molecular  scale- 
temperature  of  4400°K  Is  also  equivalent  to  any  combination  of 

m 

o 

T  »  4400°  =  —  T 
m  m 

at  this  altitude,  ranging  from  an  implausible  kinetic  temperature  of 
4400°  with  no  dissociation,  to  an  equally  implausible  314°K  in  a 
hydrogen  atmosphere.  Reality  certainly  lies  somewhere  in  between. 

The  technique's  advantage  is  this  ability  to  provide  reasonably 
correct  values  of  pressure  and  mass  density  without  precise  knowledge 
of  the  kinetic  temperature  distribution  and,  separately,  the  mean 
molecular  mass. 


TENTATIVE  STANDARD 

A  rather  high  kinetic  temperature  of  162. 5°K  was  derived  from 

(4) 

an  extension  of  the  Model  II  Mean  Atmosphere'  for  the  mesopause 
at  an  altitude  of  144  km.  Independently,  Walker  and  Jastrow  have 
recently  published  a  mesopause  temperature  for  Mars  (altitude 
unknown)  of  176°K.^^^^  From  this  level,  the  kinetic  temperature 
was  assumed  to  increase,  reaching  1000°K  at  500  km  altitude  with  a 
slight  decrease  of  mean  molecular  mass.  The  escape  level  is  reached 
between  1540  km  and  1760  km,  and  the  kinetic  temperature  was  assumed 
to  remain  constant  at  1100°K  from  1500  km  to  2800  km,  with  the 
molecular  mass  decreasing  from  16  (near  the  escape  level)  to  14 
(at  about  2800  km) . 

Any  of  these  assumptions  for  probable  limiting  conditions 
would  appear  equally  reasonable  on  the  basis  of  Chamberlain’s 
investigation  of  the  probable  mechanisms  responsible  for  the  thermal 
behavior  of  the  Martian  upper  atmosphere.  But  to  determine  the 
limiting  parametric  envelopes,  the  specific  numerical  choices  are 
expected  to  represent  also  realistic  assessments  of  the  possible 
variability.  This  will  also  roughly  account  for  the  expectation, 
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/o) 

based  on  knowledge  of  the  Earth's  upper  atmosphere^'  '  that  the 
night  temperature  at  the  thermopause  and  in  the  thermosphere  will 
be  substantially  lower  than  day  temperatures. 
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Appendlx  E 

CONSTRUCTION  PARAMETERS 


Figures  1,  2,  3,  and  Table  1  show  the  final  numerical  results  as 
the  three  tentative  model  atmospheres.  Actual  conditions  at  any 
altitude  level  should  fall  between  the  tabulated  value  given  by  the 
Tentative  Maximum  and  that  given  by  the  Tentative  Minimtim.  It  must 
be  emphasized,  however,  that  the  atmospheric  parameters  are 
tabulated  to  three  or  four  significant  figures  for  mathematical 
convenience  only.  The  spread  in  values  between  Maximum  and  Minimum 
indicates  the  uncertainty  at  any  altitude  and,  if  wished,  can  be 
interpreted  as  the  probable  error  of  the  Tentative  Standard  Atmos¬ 
phere  . 

Table  5  lists  the  principal  construction  parameters  used  in 
the  computational  program,  and  some  pertinent  data  for  different 
altitude  levels  of  special  interest.  The  reason  for  the  choice  of 
certain  input  parameters  pertaining  to  the  lower  atmosphere,  and 
planetary  parameters,  has  been  explicitly  discussed  in  an  earlier 
publication.'  ' 

Table  6  provides  information  about  the  thermal  structure  of  the 
model  atmosphere  in  terms  of  molecular  scale  temperature  and 
potential  height.  For  illustrative  purposes,  it  also  shows  one 
possible  choice  of  the  combined  variation  of  kinetic  temperature  and 
molecular  mass  for  the  given  molecular  scale  structure. 


Table 


Number  density  cm  1.487  x  10  2.463  x  10  3.20 


Table  5  (Continued) 
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Table  6 

THERMAL  STRUCTURE  OF  MODEL  ATMOSPHERES 


Gradient 

Molecular 

of  Molecular 

Representa- 

Representa- 

Potential 

Scale  Tern- 

Scale  Tem- 

tive  Kinetic 

tive  Mean 

Height 

perature 

perature 

Temperature 

Molecular 

(km) 

(°K) 

(°K/km) 

(°K) 

Mass 

TENTATIVE  MINIMUM 

0 

200 

-3.79 

200 

30 

26 

101.46 

0 

101.5 

30 

89 

101.46 

101.5 

30 

+2.315 

436.655 

906 

+1.140 

906 

30 

1044.534 

1599 

A 

906 

17 

1500. 0 

1599 

U 

906 

17 

TENTATIVE  STANDARD 

0 

250 

-3.75 

250 

29 

10 

212.5 

212.5 

29 

-2.5 

30 

162.5 

162.5 

29 

144 

162.5 

u 

+  2.99 

162.5 

29 

436 

1035.714 

1000.0 

28 

+  1.59 

1038 

1993.749 

1100.00 

16 

+0.62 

1500 

2278.569 

1100. 00 

14 

TENTATIVE  MAXIMUM 
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Appendlx  F 
EXOSPHERE 


The  validity  of  the  equations^  collected  in  Appendix  k,  becomes 

extremely  doubtful  in  the  higher  regions  of  the  exosphere.  Without 

a  magnetic  fields  the  radial  extent  of  any  planetary  exosphere^ 

rotating  with  the  planet,  has  for  its  absolute  outer  limit  the 

distance  at  which  gravitational  and  centrifugal  forces  balance. 

This  occurs  at  about  35,800  km  above  the  Earth's  equator,  and  at 

about  17,000  km  above  Mars'.  Long  before  this,  however,  the 

particle  density  of  the  Martian  exosphere  will  have  decreased  to 

that  of  interplanetary  space.  This  is  evident  from  studies  of  the 

(26) 

probable  exospheric  regimes  on  Earth  published  by  Grimminger,  '  ' 

Milne, Opik  and  Singer, Shen, and  others. 

Nevertheless,  the  previously  discussed  computer  equations 

permit  the  extention  of  numerical  integration,  with  assumptions 

about  the  radial  distribution  of  mean  molecular  mass,  up  to 

altitudes  in  excess  of  one  planetary  diameter.  While  no  great 

reliability  can  be  assigned  to  the  numerical  results,  exospheric 

density  data  thus  obtained  appear  to  fall  between  values  derived  by 
(29) 

Shen  for  a  collisionless  exosphere,  and  a  barometric  formula 
distribution,  with  constant  kinetic  temperature  or  molecular  mass, 
respectively.  The  results  indicate  the  relative  importance  of  absorp¬ 
tion  of  solar  radiation  by  different  molecular  and  atomic  species, 
even  though  they  are  quantitatively  inconclusive  as  regards  the 
actual  radial  extent  of  the  Martian  exosphere. 

For  a  special  case,  hcwever,  of  the  distribution  of  atmospheric 
mass  density  in  the  equatorial  plane,  the  calculated  upper  and  lower 
possible  extreme  limits  indicate  that  Phobos  may  conceivably  be 
subject  to  some  extremely  small  drag  effects.  Work  on  this  problem 
is  in  process. 
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